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We revise and extend the stochastic approach to cumulative weak lensing (hereafter the sGL 
method) first introduced in Ref. [1]. Here we include a realistic halo mass function and density 
profiles to model the distribution of mass between and within galaxies, galaxy groups and galaxy 
clusters. We also introduce a modeling of the filamentary large-scale structures and a method 
to embed halos into these structures. We show that the sGL method naturally reproduces the 
weak lensing results for the Millennium Simulation. The strength of the sGL method is that a 
numerical code based on it can compute the lensing probability distribution function for a given 
inhomogeneous model universe in a few seconds. This makes it a useful tool to study how lensing 
depends on cosmological parameters and its impact on observations. The method can also be used 
to simulate the effect of a wide array of systematic biases on the observable PDF. As an example we 
show how simple selection effects may reduce the variance of observed PDF, which could possibly 
mask opposite effects from very large scale structures. We also show how a JDEM-like survey could 
constrain the lensing PDF relative to a given cosmological model. The updated turboGL code is 
available at turboGL.org. 

PACS numbers: 98.62. Sb, 98.65.Dx, 98.80.Es 



I. INTRODUCTION 

Inhomogeneities in the large-scale matter distribution 
can in many ways affect the light signals coming from 
very distant objects. These effects need to be understood 
well if we want to map the expansion history and deter- 
mine the composition of the universe to a high precision 
from cosmolgical observations. In particular the evidence 
for dark energy in the current cosmological concordance 
model is heavily based on the analysis of the apparent 
magnitudes of distant type la supernovae (SNe) [2, 3]. 
Inhomogeneities can affect the observed SNe magnitude- 
rcdshift relation for example through gravitational lens- 
ing, in a way which essentially depends on size and com- 
position of the structures through which light passes on 
its way from source to observer. 

The fundamental quantity describing this statistical 
magnification is the lensing probability distribution func- 
tion (PDF). It is not currently possible to extract the 
lensing PDF from the observational data and we have to 
resort to theoretical models. Two possible alternatives 
have been followed in the literature. A first approach 
(e.g. Ref. [4-7]) relates a "universal" form of the lens- 
ing PDF to the variance of the convergence, which in 
turn is fixed by the amplitude of the power spectrum, 
cr$. Moreover the coefficients of the proposed PDF are 
trained on some specific N-body simulations. A second 
approach (e.g. Ref. [8, 9]) is to build ab-initio a model 
for the inhomogeneous universe and directly compute the 
relative lensing PDF, usually through time-consuming 
ray-tracing techniques. The flexibility of this method is 
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therefore penalized by the increased computational time. 

In Ref. [I] we introduced a stochastic approach to 
cumulative weak lensing (hereafter sGL method) which 
combines the flexibility in modeling with a fast perfor- 
mance in obtaining the lensing PDF. The speed gain 
is actually a sine-qua-non for likelihood approaches, in 
which one needs to scan many thousands different mod- 
els (see Ref. [10]). The sGL method is based on the weak 
lensing approximation and generating stochastic config- 
urations of inhomogeneities along the line of sight. The 
major improvements introduced here are the use of a real- 
istic halo mass function to determine the halo mass spec- 
trum and the incorporation of large-scale structures in 
the form of filaments. The improved modeling together 
with the flexibility to include a wide array of system- 
atic biases and selection effects makes the sGL method a 
powerful and comprehensive tool to study the impact of 
lensing on observations. 

We show in particular that the sGL method, endowed 
with the new array of inhomogeneities, naturally and ac- 
curately reproduces the lensing PDF of the Millenium 
Simulation [11, 12]. We also study a simple selection ef- 
fect model and show that selection biases can reduce the 
variance of the observable PDF. Such reduction could 
at least partly cancel the opposite effect coming from 
large scale inhomogeneities, masking their effect on the 
observable PDF. We also show how a JDEM-like sur- 
vey could constrain the lensing PDF relative to a given 
cosmological model. Along with this paper, we release 
an updated version of the turboGL package, which is a 
simple and very fast Mathematica implementation of the 
sGL method [13]. 

This paper is organized as follows. In Section II we 
introduce the cosmological background, the generic lay- 
out of inhomogeneities and review the basic formalism 
needed to compute the weak lensing convergence. In Sec- 
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tion III we derive the halo mass function and the halo 
density profiles and define the precise modeling of fila- 
ments. In Section IV we present the revised and extended 
sGL method. The exact discretization of the model pa- 
rameters, which is a crucial step in the sGL model build- 
ing, is explained in Section V and in Section VI we ex- 
plain how the realistic structures where halos are confined 
in filaments are modelled in the sGL method. Section VII 
contains our numerical results including the comparison 
with the cosmology of the Millennium Simulation [14] 
and, finally, in Section VIII we will give our conclusions. 



II. SETUP 

A. Cosmological Background 

We consider homogeneous and isotropic Fricdmann- 
Lemaitre-Robertson- Walker (FLRW) background solu- 
tions to Einstein's equations, whose metric can be written 
as: 

ds 2 = -c 2 dt 2 + a 2 {t) [dr 2 + f 2 K {r)dtt 2 } , (1) 

where dVi 2 = d9 2 + sin 2 8d(f> 2 and 

f X- 1 /2 sin ( if i/2 r ) K>0 

Mr) ={ r K = , (2) 

[ {-K)- 1 ' 2 sinh[(-Jf )V2 r ] K < 

where K/a 2 (t) is the spatial curvature of any t— slice. 

In particular we will focus on wCDM models whose 
Hubble expansion rate depends on redshift according to: 



^ ee E 2 ( Z ) = n Q0 (l + z)^ (3) 

+ n K0 (i + zf + n M0 (l + zf + n R0 (1 + zf , 

where q(z) is given by: 

q(z) = — — r / — —dz . 

qy ' ln(l + z) 7 1 + z' 



(4) 



Here w(z) could be taken to follow e.g. the parameteri- 
zation [15]: 



w(z) =W +W a j 



for which: 



q(z) = 1 + w + w a 



+ z 



w a z 



(l + z)ln(l + z) ' 



(5) 



(6) 



For a constant equation of state w(z) — w , the latter 
reduces to q(z) = 1 + wq. 

^Q0j ^mo and Qro are the present-day density param- 
eters of dark energy, matter and radiation and Qko = 
1 — £1qo — ^mo — Qrq represents the spatial-curvature 
contribution to the Fricdmann equation. Introducing the 



Hubble radius Lh = c/H and the spatial-curvature ra- 
dius Lk = a/(±K) 1,>2 , we have at any time the relation 
fix = —L 2 H /L 2 K . We will also need the matter and dark 
energy density parameters at a given time or redshift: 



(1 + z) 3 
E 2 (z) ' 
(1 + z) 3 ^ 



(7) 
(8) 



Throughout this paper, the subscript will denote the 
present-day values of the quantities. For example the 
critical density today is p C o = 3H 2 /8irG where H Q = 
lOOh km s" 1 Mpc"\ while p c = 3H 2 /8nG is the critical 
density at any time. 

Substituting in Eq. (3) H = a{t)/a{t) and 1 + z = 
ao/a(t) we obtain the equation we have to solve in or- 
der to find the time evolution of the scale factor a(t), 
the only dynamical variable in an FLRW model. We fix 
the radiation density parameter to = 4.2 • 10~ 5 /i~ 2 . 
Moreover Q,m = £Idm + ^b, that is, dark and baryonic 
matter contribute together to the total matter density. 
The line-of-sight and transverse comoving distances are: 



d C \\(z) ee r{z) = — ( 
«o Jo 

dc±(z) = f K (r(z)), 



dz' 



(9) 
(10) 



from which we find the angular and luminosity distances 
together with the distance modulus: 



d A (z) = aoil + z^fKiriz)) 
d L {z) = a (l + z)f K (r(z)) 
d L (z) 



m(z) = 5 log 



10 



10pc 



(11) 
(12) 

(13) 



For scales smaller than the spatial-curvature radius 
(±K) 1 / 2 Ar = aAr/L K < 1, f K (Ar) ~ Ar and we can 
use the Euclidean geometry. 



B. Matter Inhomogeneities 

The aim of this paper is to compute statistical weak 
lensing corrections to the measured light intensities, gen- 
erated by the inhomogeneous matter distribution along 
the line of sight to a distant object. The basic quantity 
we need to model is the matter density contrast 5m- 



&M{r,t) = 



Pm(r,t) 

PM(t) 



- 1 



(14) 



where the lowercase p m indicates the local and inhomo- 
geneous matter field while pm — &m Pc is the time de- 
pendent average mass density. The density contrast di- 
rectly enters the expression for the weak lensing conver- 
gence [16]: 

k= / drG(r,r s ) pmc Sm , (15) 
Jo 
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where r s is the co-moving position of the source and the 
integral is along an unperturbed light geodesic. The den- 
sity Pmc = a o Pmo is the constant matter density in a 
co-moving volume and we defined the auxiliary function 1 

G(r ,, ) = ^M2£gpi>i. (16) 

which gives the optical weight of an inhomogeneity at 
the comoving radius r. The convergence is related to the 
shift in the distance modulus by: 

Am = 51og loA1 - 1 / 2 ~51og 10 (l-«), (17) 

where p is the net magnification and the second-order 
contribution of the shear has been neglected [1, 16]. It 
is obvious that an accurate statistical modeling of the 
magnification PDF calls for a detailed description of the 
inhomogeneous mass distribution. 

In this paper we will significantly improve the modeling 
of the inhomogeneities from our previous work [1] , where 
only single-mass spherical overdensities were considered. 
First of all, we improve the modeling of these "halos" 
by using a realistic halo mass function f(M,z), which 
gives the fraction of the total mass in halos of mass M 
at the redshift z. The function /(M, z) is related to the 
(comoving) number density n(M, z) by: 



dn(M,z) = n(M,z)dM 



Pmc 

"AT 



f(M,z)dM, (18) 



where we defined dn as the number density of halos in 
the mass range dM. The halo function is by definition 
normalized to unity 



/ 



f(M, z)dM = 1 



(19) 



The idea is of course that halos describe large virialized 
mass concentrations such as large galaxies, galaxy clus- 
ters and superclusters. Not all matter is confined into 
virialized halos however. Moreover, only very large mass 
concentrations play a significant role in our weak lensing 
analysis: for example the stellar mass in galaxies affects 
the lensing PDF only at very large magnifications [11, 12] 
where the PDF is close to zero. Therefore the fraction of 
mass A./h concentrated in large virialized halos can be 
defined by introducing a lower limit to the integral (19): 



AM*) = f 



f{M, z)dM < 1 



(20) 



Only mass concentrations with M > M cut are treated 
as halos. The remaining mass is divided into a family 
of large mass, but low density contrast objects with a 



fraction AJl and a uniform component with a fraction 
A/t/ = Pmu/pm, such that 



Af H + Af L + Af u = l. 



(21) 



The low contrast objects are introduced to account for 
the filamentary structures observed in the large scale 
structures of the universe. In our analysis they will be 
modeled by elongated objects with random positions and 
orientations. The mass in these objects can consist of a 
smooth unvirialized (dark) matter field and/or of a fine 
"dust" of small virialized objects with M < M cut . In 
weak lensing this distinction does not matter because 
small halos act effectively as a mean field, with a size- 
able contribution only at very large magnifications (see 
comment above about stellar mass). 

For later use we define Afm = AJl + Afu = 1 — Afn 
which gives the total mass fraction not in large virialized 
halos. If we only consider virialized masses larger than 
companion galaxies (M cut ~ 10 10 /i -1 M ), then typical 
values for the concordance model are Afm ~ 0.5 at 
z = 0[14] and A/lu ~ 0.7 at z — 1.5, with a weak depen- 
dence on the particular /(M, z) used. Although Afm is 
sensitive on M cut , the lensing PDF depends only weakly 
on the cut mass. Moreover, halo functions obtained from 
N-body simulations are valid above a mass value imposed 
by the numerical resolution of the simulation itself and 
so the use of M cut is also necessary in this case [17]. 

To obtain an as accurate modeling as possible, we 
will use realistic mass functions and spatial density pro- 
files for the halo distribution. There is less theoretical 
and observational input to constrain the mass distribu- 
tion of the filamentary structures or their internal den- 
sity profiles. We will therefore parametrize the filaments 
with reasonable assumptions for their lengths and widths 
and by employing cylindrical nonuniform density pro- 
files. Our modeling allows treating the two families of 
inhomogeneities independently. Both can be given ran- 
dom spatial distributions, or alternatively all halos can 
be confined to have random positions in the interiors of 
randomly distributed cylinders. The latter configuration 
more closely resembles the observed large scale struc- 
tures, and an illustration created by a numerical sim- 
ulation using the turboGL package is shown in Fig. 1. In 
Section VI A we will discuss the power spectrum and the 
large-scale correlations of our model universe. 

We can now formally rewrite Eq. (14) for our model 
universe in which the density distribution is given 

by Pm = J2j P] + PMU- 



Sm = 



Pmc 



+ Afu-l 



(22) 



1 Note that this definition slightly differs from the one of Ref. [1]. 



where the index j labels all the inhomogeneities, we used 
Eq. (21) and we defined pj = a~ 3 Mjtpj, so that both 
virialized halos and the unvirialized objects are described 
by a generic reduced density profile tpj which satisfies the 
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FIG. 1: Shown is an illustrative-only projection of the random 
matter density within a 100/i -1 Mpc thick slice generated by 
our stochastic model. The shaded disks represent clusters and 
the shaded cylinders the filamentary dark matter structures. 
Only large clusters are displayed. 



A light ray that misses all the inhomogeneities ipj will 
experience a negative total convergence hue — k>u + k e, 
which gives the maximum possible demagnification in a 
given model universe. In an exactly homogeneous FLRW 
model the two contributions krl and kue cancel and 
there is no net lensing. In an inhomogeneous universe, 
on the other hand, a light ray encounters positive or neg- 
ative density contrasts, and its intensity will be magni- 
fied or demagnified, respectively. The essence of the sGL 
method is finding a simple statistical expression for the 
probability distribution of the quantity n HL in the inho- 
mogeneous universe described above. 

For a discussion about the validity of the weak lens- 
ing approximation within our setup see Ref. [1] , where it 
was shown that the error introduced is < 5%. We stress 
again that the lensing caused by stellar mass in galaxies 
is negligible in the weak lensing regime [11, 12] and so we 
can focus on just modeling the dark matter distribution 
in the universe. Also, we will treat the inhomogeneities 
as perturbations over the background metric of Eq. (1). 
In particular we will assume that redshifts can be related 
to comoving distances through the latter metric. See 
Ref. [18 21] for a discussion of redshift effects. 

In the next section we will give the accurate model- 
ing of the inhomogeneities. We begin by introducing the 
halo mass function and the detailed halo profiles, and 
then move on to describe the precise modeling of the 
cylindrical filaments. 

III. HALO AND FILAMENT PROPERTIES 



normalization J v fdV — l. 2 The lensing convergence we 
are interested in can now be written as 

k = k hl + k v + k e , (23) 

where the positive contribution due to the inhomo- 
geneities (halos and low contrast objects) is 

khl = f dr G{r, r„) V M jiPj , (24) 
Jo 

the also positive contribution due to the uniformly dis- 
tributed matter is 

ku = Pmc drG(r,r 3 )Afu(z(r)), (25) 
Jo 

and the negative empty beam convergence is: 

ke = ~Pmc drG(r,r s ). (26) 
Jo 



We begin by explaining our dark matter halo modeling. 
The two main concepts here are the halo mass function 
f(M, z) giving the normalized distribution of the halos 
as a function of their mass and redshift, and the dark 
matter density profile within each individual halo. Both 
quantities are essential for an accurate modeling of the 
weak lensing effects by inhomogeneities. We shall begin 
with the halo mass function introduced above in Eq. (18). 

A. Halo Mass Function 

The halo mass function acquires an approximate uni- 
versality when expressed with respect to the variance of 
the mass fluctuations on a comoving scale r at a given 
time or redshift, A(r, z). Relating the comoving scale r to 
the mass scale by M = -^r 3 pmc, we can define the vari- 
ance in a given mass scale by A 2 (M, z) = A 2 (r(M),z). 
The variance A(r, z) can be computed from the power 
spectrum: 

A\r,z)=(^ = J™ d ^A\k,z)W*{kr), (27) 



2 Note that this definition slightly differs from the one of Ref. [1]. 



where W(kr) is the Fourier transform of the chosen (top- 
hat in this work) window function and the dimensionless 
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power spectrum extrapolated using linear theory to the 
rcdshift z is: 



A 2 (M) 



2tt 2 



P(k,z) 

(—) 

\aoHoJ 



(28) 



3+n s 



T 2 (k/a )D 2 (z). 



Here n s is the spectral index and Sho is the amplitude 
of perturbations on the horizon scale today, which we fix 
by requiring: 



A(r = 8/a h Mpc, z = 0) = <t 8 . 



(29) 



where the value of cr 8 is estimated by cluster abundance 
constraints [22]. D{z) is the linear growth function which 
describes the growth speed of the linear perturbations in 
the universe. Fit functions for D could be found, e.g., 
in Ref. [23], but it can also be easily solved numerically 
from the equation: 



D"(z) + lf^[n M (z) + (3w(z) + l)n Q (z)] 
I 1 + z 



3 D(z) 
2jl + zj- 



£Im(z) = 0, 



(30) 



where we have neglected the radiation. As usual, we nor- 
malize D to unity at the present time, D(0) — 1. Finally, 
for the transfer function T(k) we use the fit provided by 
the Equations (28-31) of Ref. [24], which accurately re- 
produces the baryon-induced suppression on the interme- 
diate scales but ignores the acoustic oscillations, which 
are not relevant for us here. 

With A(M, z) given, we can now define our halo mass 
function. Several different mass functions have been in- 
troduced in the literature, but here we will consider the 
mass function given in Eq. (B3) of Ref. [17]: 

fj(M,z) = 0.301 exp(-|lnA(M,z)- 1 +0.64| 3 - 82 ), (31) 

which is valid in the range —0.5 < In A -1 < 1.0. Because 
of the change of variable, fj is related to our original 
definition of / by: 



f(M,z) = fj(M, z) 



dlnA(M,z)" 
dM 



(32) 



The mass function of Eq. (31) is defined relative to 
a spherical-overdensity (SO) halo finder, and the over- 
density used to identify a halo of mass M at redshift z 
is Ago = 180 with respect to the mean matter density 
Pm{z). The SO finder allows therefore a direct relation 
between halo mass M and the radius R v : 



M = —R 3 pPM (z)A so . 



(33) 



The subscript p will denote the proper values of otherwise 
comoving quantities throughout this paper. For example, 
the comoving halo radius R is related to the proper value 
by R p — a R. 



A direct relation between the mass and the radius of 
a cluster is necessary for our sGL method. This is why 
we prefer mass functions based on SO finders over mass 
functions based on friends-of-friends halo finders; for the 
latter an appropriate Ago to be used in Eq. (33) is not di- 
rectly available. Moreover, as shown in Ref. [17, 25], the 
SO(180) halo finder gives a good degree of universality 
to the mass function. Another mass function which man- 
ifests approximate universality with the SO (180) halo 
finder [25] is provided by Sheth & Tormen in Ref. [26]. 



B. Halo Profile 

With the halo mass function given the only missing 
ingredient is the halo profile which, as said after Eq. (22), 
we describe by means of the reduced halo profile tph- We 
stress that our halo of a given mass M and redshift z 
is an avererage representative of the total ensemble of 
halos, which in reality have some scatter in the density 
profiles. 

We will focus on the Navarro-Frenk- White (NFW) pro- 
file [27]: 



Ph 



Pm {r p /R 8 ){l + r p /R s y 



(34) 



Here the scale radius is related to the halo radius by 
R p = cR Sl where c is the concentration parameter to 
be defined below. By integrating equation (34) one finds 
the total mass M = 4ir 5 C p M (R P /c) 3 [ln(l + c) - c/(l + c)} 
which, when combined with Eq. (33), gives: 



A 



so 



3 ln(l + c) - c/(l + c) 



(35) 



This relation fixes S c once the concentration parameter 
is given. The reduced halo profile <ph in comoving coor- 
dinates is then: 



¥>h{r,t) 



, , R 

iirr I — + r 



ln(l + c) 



1 + c 



(36) 

where the ^-dependence (as well as M-dependence) arises 
through the redshift dependence of R (Eq. (33)) and c 
(Eq. (37)). 

The last ingredient needed to fully specify the NFW 
profile is a relation between the concentration parameter 
and the halo mass at a given redshift. For the ACDM 
model we will use the following fit obtained from numer- 
ical simulations [28] satisfying WMAP5 cosmology [29]: 



c{M,z) = 10.14 



M 



-0.081 



(1+*) 



-1.01 



(37) 



where the pivot mass scale is Mr = 2 • 10 h M Q . 
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Apart from the NFW profile, alternative simple density 
profiles of interest include: 



fh 
fh 



AttR 3 



= 0.97 



1 



.! exp[-r 2 /2(fl/3) 2 
(2tt) 3 / 2 (.R/3) 3 



4irr 2 R 



Uniform, (38) 
Gaussian, (39) 
SIS, (40) 



where R is again given by Eq. (33). The normalization 
factor in (39) is necessary to have the correct volume 
integral of unity and is due to the restriction of the profile 
between and R. In these cases there are no further 
parameters to be fixed. 



C. Filaments 

We will model dark matter filaments by cylindrical, 
non-uniform low density objects. Let us stress that 
throughout this chapter by a "filament" we refer only 
to the smoothly distributed dark matter component of a 
real filament; we shall later explain how virialized halos 
can be embedded within these objects to form a more 
realistic model of the true filamentary structures in the 
full sGL framework. Keeping the low contrast compo- 
nent and the halos as separate entities in our evaluation 
is convenient for the modeling and explicitly avoids dou- 
ble counting. 

Using Eqs. (18) and (21) we can obtain the follow- 
ing equation which relates the comoving filament density 
Anf and the filament mass Mf to the fraction Afm of 
mass not confined in virialized halos: 



Mf An f = pmc Pf Mlv , 



(41) 



where the quantity /?/ defines the fraction of the total 
unvirialized mass which is confined in filaments: A/l = 
/3fAfnj. Because of the ongoing halo formation A/lu 
decreases with time and we shall simply assume that /3f 
is a constant. 

We will assume that all filaments have the same mass, 
but this assumption can easily be relaxed later by intro- 
ducing a filament mass function [30] . We are also neglect- 
ing merging of smaller filaments into larger ones, whereby 
the comoving number density of filaments An / remains 
a constant. This assumption seems to be supported by 
numerical simulations [14] at least at low redshifts. From 
Eq. (41) it then follows that the filament mass Mf is not 
constant, but it gradually decreases as more and more 
halos virialize. (The total mass of the full filamentary 
structure which includes the halos remains on average 
constant.) 

Now we need to specify the precise geometry, dimen- 
sions and time-dependent density profile <pf for our fil- 
aments. First, we will relate the filament mass to its 
length and radius at present day according to: 



where R p o and L p o are the present-day filament proper 
radius and length and A^o is the average overdensity of 
the filament with respect to the FLRW matter density. 
Consistent with our assumptions of neglecting filament 
mergers and using a constant comoving filament density, 
we will assume that filaments have a constant comoving 
length. This is again in accordance with the generic pic- 
ture of a web of filaments of constant comoving length 
condensing with time seen in numerical simulations. Fi- 
nally, for simplicity, we will choose our filaments to have 
a constant proper radius. Summarizing, we will use: 



L 

R(t) 



Lpo/ao , 
R p0 /a(t) . 



(43) 



The simplest assumption would be to take filaments to 
have uniform density. In this case the reduced density 
profile is simply iff = Vf~ . However, we will also con- 
sider a profile which is uniform along the length of the 
filament L, but which has a gaussian profile in the radial 
direction. In this case 



ip f (r,t) = 0.989" 



exp[-r 2 /2(fl/3) 2 ] 
L2tt(R/3) 2 



(44) 



where < r < R and < I < L, and tpf = otherwise. 

Finally, let us compute the mass of the "dressed" fil- 
aments which comprises the "bare" filaments until now 
discussed together with the halos (to be embedded in 
Section VI). It is easy to find that: 



Mf Q 



Vf P o Pmo 



1 



A/l 



(45) 



where also the contribution from the uniformly dis- 
tributed matter is included. Eq. (45) is valid if all halos 
are embedded within the filaments and has to be slightly 
modified if different levels of confinement are considered. 

This completes our description of the objects in our 
model. Next we shall have to understand how to simulate 
their effect on the weak lensing properties of the complete 
inhomogeneous model universe. 



IV. THE sGL METHOD 

In this section we will review and extend the stochastic 
gravitational lensing (sGL) method of Ref. [1] for com- 
puting the probability distribution function of the lens 
convergence k in the presence of inhomogeneities. The 
basic quantity to evaluate is the inhomogeneity-induced 
part Khl in Eq. (24): 



khl(zs) 



3 



rjlW, (46) 



M 



fo = 7ri? 2 L p0 Pmo A/ = Vf p0 Pmo A/o , (42) 



where r s = r(z s ). We wish to obtain a probabilistic pre- 
diction for this quantity along a line of sight to a source 
located at r s , through a random distribution of halos and 
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filaments. The problem could be solved by constructing a 
large comoving volume of statistically distributed objects 
such as the one shown in Fig. 1, and then computing the 
distribution of khl along random directions in this space. 
However, an equivalent and computationally much more 
efficient approach is to construct random realizations of 
object locations along a fixed geodesic and compute the 
lensing PDF from a large sample of such realizations. 

Let 9 refer to a particular realization of the integral in 
Eq. (46) and denote the resulting convergence by k(9, z s ). 
Because of the finite size of the halos and filaments, only 
a finite number Ng of objects intercept the geodesic and 
contribute to the sum. Moreover, because all objects are 
small compared to a typical comoving distance to the 
source r s 3> R, L, the function G is, to a good approxima- 
tion, a constant G(r,r s ) w G(rj,r s ) for each individual 
object. Similarly, one can assume that ip(x,t) w tp(x,tj) 
with tj — t(rj), and so one finds: 

N a 

KH^z^^J^Girj^^itj), (47) 

where Sj is the surface mass density of the j-th in- 
tercepting object in the realization. Next divide the 
geodesic into N$ subintervals with centers at f, and 
with (possibly variable) length Ar^ -C r s , such that 
G(r,r s ) « G(n,r s ) = G,(r s ) = Ar' 1 G(r, r s )dr 
holds within each interval. In this way we can write 

N s N f 

K HL {6,z s )^J^J^G i {r s )^ u {t i ), (48) 
»=i jt=i 

where fi labels all objects encountered by the geodesic 
within the comoving length bin An, such that J2fJi N% — 
Ng . At this level Eq. (48) still consists of a sum over all 
objects along a given path. Next we categorize these ob- 
jects into different classes depending on the parameters 
that define the surface density Ej 4 . For spherical halos 
of a given (^-profile X depends only on the mass of the 
halo and the impact parameter, and for filaments the 
relevant parameters are the angle between the main fil- 
ament axis and the geodesic and the impact parameter 
in the cylindrical radius. Now assume all these variables 
are discretized into finite-length bins. 3 The precise way 
of the discretization will be discussed in detail in Section 
V, but for now we simply let a generic index u label the 
set of independent parameter cells created by the bin- 
ning. If we now let k iu denote the number of objects 
that fall into the cell labeled by indices i and u in the 
realization 9, we can replace the sum over the individual 



objects in Eq. (48) by a weighted sum over the discrete 
parameter cells: 

N s N c 

k hl (9,z s ) ~ Y^Giir^Y^hlXiuiU) 

i—l u—1 

N S N C 

i— 1 u—1 

where Nq is the total number of independent parameter 
cells 4 used in the binning, Y^fJi l k L = N e and K iiu 
is the convergence due to one object in the bin iu: 

Ki iu (z s ) = Gi(r s )T, iu (ti). (50) 

Note that the quantity Ku u is a function of the distance 
and internal variables, it is universal for arbitrary real- 
izations: all the information specific to a particular real- 
ization 9 is contained in the set of integers {kf u }. 

Equation (49) is the starting point of the sGL analy- 
sis, because it can be easily turned into a probabilistic 
quantity. Indeed, instead of thinking of realizations 9 
along arbitrary lines of sight through a pre-created model 
universe, we can define a statistical distribution of con- 
vergences through Eq. (49). Indeed, we have shown in 
Rcf. [1] that, for random initial and final points of the 
geodesies, the integers kf u are distributed as Poisson ran- 
dom variables: 

P kiu = {A ^ U p 1 e-^, (51) 

where the parameter AiVj M is the expected number of ob- 
jects in the bin volume AVi U : 

AN iu = An lu AV iu = An m An AA iu . (52) 

Here Ani U is the comoving density of objects correspond- 
ing to the parameters in the bin iu and AA iu is the cor- 
responding cross sectional area of these objects in co- 
moving units. We shall specify these quantities precisely 
in Section V. Physically, the statistical distribution of 
convergences is equivalent to our original set of realiza- 
tions 9 averaged over the position of the observer. This 
is a welcome feature because the statistical model explic- 
itly incorporates the Copcrnican Principle. The original 
realization 9 has now been replaced by a configuration 
of random integers {k iu } and the convergence equation 
Eq. (49) by a statistical convergence: 

K HL{{kiu}, Z s ) = ^ ' k; lu K\i u 

(z s ) , (53) 

iu 



3 All other parameters, such as filament length, radius and mass 
are kept fixed. In a more accurate modeling the fixed parameters 
(for example the filament mass) could be let vary, in which case 
they should also be discretized. 



4 The binning of internal variables could depend on r, and so the 
indexing u and the number of different parameter cells Nq could 
in fact depend on i, but we suppress this notation for simplicity. 
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where the probability for a particular configuration to 
occur is just 

N s N c 

p {^} = n n p ^ • ( m ) 

i— 1 u—1 

It is easy to show that this probability distribution is cor- 
rectly normalized. Moreover, including the convergence 
kjje due to the uniformly distributed matter and the 
empty space (see Eq. (26) and below) , one can show that 
the total convergence corresponding to the configuration 
{ki U } becomes [1]: 

K({k iu }, z 8 ) = ^2 K nu(z s ) (kiu - AN iu ^j . (55) 

iu 

For a configuration without any inhomogeneities we cor- 
rectly recover n{{k iu — 0},z s ) — kjje{z s ), see Section 
V. Moreover, Eq. (55) shows explicitly that the expected 
convergence vanishes consistently with the photon con- 
servation in weak lensing, because (k iu ) — AN iu for a 
random variable following Eq. (51). The final conver- 
gence PDF can now be formally written as 

1 /•*+# 

P w \(k, z s ) =\im A ^ — ^ dn' P w i(k' , z 8 ) , (56) 
where 

P w1 (k,z s )= ^ P {kin} S(K- K({k iu },z s )) , (57) 

is a discrete probability distribution. Note that the most 
likely configuration which maximises P{k iu } corresponds 
to the mode of the Poisson distribution, which is the 
floor of its parameter: fcj„ — > [AN iu \. Moreover, for 
a large ANi U the Poisson distribution approximates a 
gaussian with mean and variance equal to AN iu and the 
most likely configuration therefore approaches the mean: 
ki U — > ANi U . When this is the case the mode of the 
lensing PDF vanishes even for a single observation and 
the PDF tends to a gaussian. One can similarly create 
the PDF in the shift in the distance modulus Am or 
in the magnification fi, using Eq. (17) to compute these 
quantities for a given configuration 5 and then replacing 
k by the desired quantity in equations (56-57). 

In practice, equations (56-57) arc not useful for a direct 
evaluation of the PDF. Indeed, the configuration spaces 
are infinite and a direct evaluation of the sum of configu- 
rations would not be feasible even for relatively roughly 
discretized systems. Instead, we compute an approxima- 
tion for P w i by statistically creating a large set of ran- 
dom configurations {k iu } from the probability distribu- 
tion P{k iu }, and by forming a discrete normalized his- 
togram out from the corresponding set of convergences 



(or distance moduli or magnifications). That is, we define 
a histogram PDF using the function: 

P wl (Ri,z a )Ai»^^, (58) 

J *sim 

where iV sun is the total number of realizations and iV s i mji 
is the number of convergences in the sample falling in 
the bin A; centered around the mean value Ri. It is 
easy to check that Eq. (58) is correctly normalized to 
unity. Moreover, the more configurations one creates in 
the simulation step, the more fine-detailed and accurate 
approximation for P w \(k, z s ) can be constructed. 6 

To summarize, the sGL method consists of two steps: 
first the internal variables specifying the inhomogeneities 
are binned, reducing them to a finite number of different 
object classes labeled by the cell index u, for which uni- 
versal convergence functions n liu can be computed. Sec- 
ond the occupation numbers ki U of these objects within a 
given co-moving distance bin are shown to be Poisson dis- 
tributed with parameter ANi U . The lensing PDF is then 
computed statistically, binning the convergence distribu- 
tion obtained from a large set of random configurations 
{k iu } on a given geodesic. Computing the lensing PDF 
this way is very efficient. Finding the PDF for a given 
model universe using the turboGL package [13] typically 
takes a few seconds in an ordinary desktop computer. See 
the Appendix A for more details about its performance. 

Let us point out that Eq. (58) treats all objects on the 
same footing. That is, all object families are indepen- 
dently randomly distributed throughout the space. Con- 
fining the halos to within filaments will require slight 
modifications which we will explain in the section VI. 
For now, we shall study the effects that may cause the 
actual observed PDF deviate from the fundamental weak 
lensing PDF. 

A. Observable PDF 

The lensing PDF (58) is the fundamental quantity in- 
herent to a given background model and a set of inho- 
mogeneities. However, in reality there are interferences 
that prevent us from observing P w \ directly. Firstly, the 
intrinsic magnitudes of the sources are not known with 
arbitrary accuracy, and so the observed lensing PDF is 
at best a convolution: 

P 1 (Am,z s )= J dyP wl (y,z s )P in (Am-y), (59) 

where P; n describes the source magnitude dispersion for 
an imperfect standard candle. Clearly the fundamental 



This is true assuming that nonlinear effects can be neglected, 
i.e., in the weak lensing limit. 



One could also define a smooth P w i(k, z s ) from the simulated 
data as a moving average weighted by some window function, but 
that would not give any advantage in the analysis with respect 
to using the simple histogram of Eq. (58). 
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lensing PDF is recovered in the limit that P; n becomes a 
delta function. If there were no other sources of errors, 
the distribution p could be used to contrast observations 
against different cosmological background models, sets of 
inhomogcneities and the models for the source magni- 
tude distributions. However, there can be several differ- 
ent types of selection effects caused by matter between 
the source and the observer or by the search strategies, 
that may bias the observed PDF from the one predicted 
by Eq. (59). One of the strengths of the sGL method is 
that it can be easily extended to model the effects of a 
wide range of selection biases on the observable PDF. 



Observational Biases 

Let us begin with a simple class of biases that might 
affect the observed distance modulus beyond the weak 
lensing correction; an obvious example would be an error 
made in the estimate of reddening. What makes such 
correction nontrivial, is that one would expect the error 
to be correlated with the types of mass concentrations 
encountered by the photon beam on its way from the 
source to the observer. The effects should thus be differ- 
ent along different lines of sight, affecting mostly those 
directions where light has to travel through volumes with 
higher mass densities. In the sGL method such correlated 
biases are very easily modelled, as one needs merely to 
replace 

Am -> Am({k iu },z s ) + 8m({k iu },z s ) 

= Am cS ({k iu },z s ) , (60) 

where the first part Am({k iu } 7 z„) is the usual weak lens- 
ing correction and the second part the object-specific cor- 
rection. Indeed, the bias function Sm({ki U }, z s ) can be 
correlated with the redshifts z s and Zi, and with all the 
internal characteristics of the objects, such as the halo 
and filament masses and impact parameters. A quantita- 
tive analysis of the possible magnitude of the function 5m 
would require detailed astrophysical information and it is 
beyond the scope of this paper. However, once properly 
modelled, the sGL allows a very easy way of evaluating 
the outcome of biases on the observations. 

Other types of selection biases might not directly af- 
fect the magnitudes, but rather the relative likelihood 
of observing particular events, in a way that is corre- 
lated with the mass density. This includes, for exam- 
ple, all sources leading to obscuration of the light beam, 
either alone or in combination with restrictions coming 
from imperfect search efficiencies and strategies. For ex- 
ample selection by extinction effects or by outlier rejec- 
tion mainly relate to high-magnification events which are 
clearly correlated with having high intervening mass con- 
centrations (halos with large M and/or small b) along the 
light geodesic. Similarly, short duration events might be 
missed by search telescopes or be rejected from the data 
due to poor quality of the light curve, for example in 



cases where a supernova is not separable from the im- 
age of a bright foreground galaxy. Probability of such 
events would obviously be correlated with the brightness 
of the source (z s and source dependent parameters), with 
the density and redshift of the intervening matter and of 
course with the search efficiencies. 

The selection biases can be studied in the sGL frame- 
work by introducing a new, object-dependent survival 
probability function. That is, we replace 

p {^} P {L} = K P ^} p {^} (61) 

in the probability distribution for the magnification 
PDF (57). Here K is a normalization constant that 
makes sure that the final PDF is normalized to unity. 
The most generic form of the survival probability func- 
tion is 

piki} = Il(p?: r ) kiu ■ (62) 

iu 

This allows the survival function to depend on the ar- 
bitrary local properties along the photon geodesic. The 
physical interpretation of P™ r is of course clear: it gives 
the relative probability that an event whose signal inter- 
cepts an object with characteristics given by internal pa- 
rameters u at redshift Zi would make it to the accepted 
observational sample. Given the form (62) the proper 
normalization is easily seen to be: 

"'ZU- 
ZU 

where 

AJVf = P™AN iu . (64) 

The correct expression for the convergence is still given 
by Eq. (55), with the important difference that the ran- 
dom integers k iu used to evaluate Eq. (58) are now drawn 
from the Poisson distribution of Eq. (63) with parameter 
the effective expected number of objects ANf^ . 

As mentioned above, modeling P|^ r quantitatively 
would need detailed input of the astrophysical proper- 
ties of the intervening matter distributions and of the 
observational apparatus, which is beyond the scope of 
this work. See Ref [1] for an analytical example giv- 
ing (k(z s )) = (I — a) k e , where a is the effective filling 
factor of a partially-filled beam [18] (see [31] for an ex- 
tension of the partially-filled beam formalism). We will 
discuss an illustrative numerical example in Section VII 
where we will consider as survival probability a simple 
step-function in the impact parameter. 

To summarise, let us note that it is easy to combine 
both types of observational biases discussed in this sec- 
tion. The true observational PDF can be defined as 

P obs (Am, z s )= Jdy P$(y, z s ) P in (Am - y) , (65) 

where p n again is the intrinsic source magnitude dis- 
persion and the effective lensing probability distribution 
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function P^(y, z s ) is computed from Eq. (58) using the 
effective probability distribution of Eq. (63), and com- 
puting the effective magnitudes for configurations from 
equations (60) and (17). This setting clearly allows mod- 
eling a very complex structure of observational biases. 

Finally note that the generic trend among all the ef- 
fects discussed above is to suppress events with higher 
magnifications. It could then be that a non-negligible 
overestimate of the observed magnitudes follows from ne- 
glecting the selection biases, even if the individual effects 
were relatively small. Moreover, while the present data 
cannot yet constrain the lensing PDF itself, it can put 
bounds on its variance. Since the observational biases are 
suppressing the high magnification tail, they give P bs a 
smaller effective variance which may open the possibility 
for a comparison with the experiments. We will further 
discuss these topics in Section VII. 



B. Binned Data Samples 

We conclude this section by deriving an effective distri- 
bution Pm q for a data sample {No}, where the set {No} 
refers to a given binning of the original data. P^ could 
be computed directly from the fundamental PDF [1], but 
it is much faster to create it directly during the initial 
simulation, bypassing the calculation of the fundamen- 
tal P w i altogether. Indeed, if we label the configurations 
within a given sample of observations by s, the mean 
convergence after No observations is: 



KNo({ k iu}) = 



1 N 

■jy— ^ n{{kiu}s) 



(66) 




where in the second line we have used the fact that the 
quantities Ku u are independent of s and in the last line 
we have used the property that the independent Poisson 
variables fcj„ iS (with the same weight) sum exactly (in 
a statistical sense) into the Poisson variable ki U ^ 01 
parameter given by the sum of the individual parameters, 
which is No AiVj„. P^ a is then given by an expression 
similar to Eqs. (56-57) or Eq. (58). 

Note, however, that including the selection effects (see 
Section IV A) in general does not commute with taking 
an average over the observations. In other words, if the 
No measurements are correlated, we cannot use Eq. (66), 
but we have to start from the fundamental PDF. On the 
other hand, Eq. (66) displays explicitly the effect of the 
size of the data sample: even if k had a skewed PDF and 
a nonzero mode for No = 1, for large No the distribution 
approaches a gaussian and eventually converges to a 5- 
function with a zero convergence, as it is clear from the 
properties of Poisson variables. 



C. Integral Formulas for Expected Value and 
Variance of the Lensing Convergence 

From the general expression Eq. (66), including a non- 
trivial survival probability, it follows that the expected 
value and variance of the convergence are: 



5> liu (P™-l)AN h 



and 



N ^ 

IU 



2 psur a at 
K liu r iu ^ Jv r< 



(67) 



(68) 



where we assumed that the survival probability for a light 
ray going through the uniformly distributed matter den- 
sity pmu equals unity and that the random variables k iu 
are uncorrelated. Eqs. (67-68) can then be put back in 
integral form giving the following exact results: 

(k) = / drG(r,r s ) dn(M,z(r))x (69) 
Jo «„, 

x / dA(b,M)(P SUT -l)X(b,M,z(r))), 
Jo 

and 

a 2 K = — / drG 2 (r,r s ) dn(M,z(r))x (70) 

N Jo JM cut 
rR(M,z(r)) 

x / dA(b,M)P sm Y?(b,M,z{r))), 
Jo 

where the integral limits for the last two integrals are 
implicitly defined and the explicit form of surface area 
element dA depends on the geometrical properties of the 
object of mass M. For example, for spherical halos it is 
dA = 2irbdb. We again stress that the quantity P sur = 
P SUI (b, M, z, z s ) is a generic function that describes the 
probability that a light ray is observed when it hits a halo 
of mass M at the redshift z and impact parameter b for 
a source at redshift z s . 

Eqs. (69-70) are a simple and direct prediction of the 
sGL method and allow to draw some general considera- 
tions. First, if P sur = 1, the expected value of the conver- 
gence is correctly zero, showing the "benevolent" nature 
of weak lensing corrections for unbiased observations. If 
however the survival probability is not trivial, selection 
biases persist even in very large datasets so that the av- 
erage convergence over a large number of observations 
approaches a nonvanishing value. Second, Eq. (70) is a 
product of positive quantities and so a nontrivial survival 
probability always reduces the observed variance. 

We would like to stress that it is easy to numerically 
evaluate the integrals of Eqs. (69-70), and their predic- 
tions can be straightforwardly implemented in \ 2 analy- 
ses based on gaussian likelihoods, which may be a reason- 
able approximation if the observable PDF is not strongly 
skewed. 
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V. DISCRETIZING THE MODEL 
PARAMETERS 

The key role played by the discretization of the model 
parameters in the sGL method is now quite obvious. Of 
course, the final lensing PDF should not depend on the 
discretization. Indeed, as the numerical value of a Rie- 
mann integral does not depend on the binning chosen if 
the integrand does not significantly vary within the bins, 
the resulting PDF will not depend on the specific binning 
if the cosmological functions Ku u are approximately con- 
stant within the bin cells iu. Once this condition is well 
met a further refinement of binning does not increase the 
accuracy. This statement can be proved formally as fol- 
lows: imagine that we have a valid binning prescription 
A and we want to move to a new prescription B which 
is N times finer. Because A is valid, the convergences 
K iiu are approximately constant within each ,4-bin. The 
relative Poisson variables can then be resummed as done 
in Eq. (66), and we obtain the same formal expression in 
terms of ,4-bins we started from. 

This theorem also explains why we do not need to bin 
any variable in which the surface mass density E of a 
given object is constant, such as the rotational angle 
along the center of a spherical halo in the lens plane. As 
the convergences Ku u are constants within angular bins, 
their contributions can be resummed to give a model 
without any angular binning. That is, the fundamental 
halo objects are actually finite- width rings in the conver- 
gence plane. 

The first quantity to be binned is the co-moving dis- 
tance to the source. Here the basic quantity that needs 
to be accurately modelled is the function G(r,r s ), which 
gives the optical weight of a given convergence plane. 
This is a rather smooth function with a shallow peak, 
and only slightly dependent on r s when expressed in the 
scaled variable r/r s . Simple linear bins fj have been 
found to give accurate results, with only 10 — 15 redshift 
slices. This turns out to be the typical number of bins 
needed for an accurate modeling of most of the internal 
variables as well. 



A. Discretizing the Halo Parameters 

We begin by discretizing the halo parameters, which 
are the mass and the impact parameter in the lens plane. 
Binning the halo masses involves two main considera- 
tions. First, the halo mass function only carries informa- 
tion about high-contrast virialized halos within a specific 
mass range (see Eq. (31)). Moreover, as explained in Sec- 
tion II B, only very large mass concentrations play a sig- 
nificant role in our weak lensing analysis. Consequently, 
we introduced the lower cut M cut to the mass distribu- 
tion so that only mass concentrations with M > M cut are 
treated as halos, while the remaining mass has been di- 
vided into low contrast objects and uniform matter den- 
sity (see Eq. (21)). Second, the fact that the concentra- 




FIG. 2: Shown is a comoving segment of a photon geodesic 
intercepting a number of halos represented by shaded discs. 



tion parameter of Eq. (37) follows a power law suggests 
the use of a logarithmic binning in mass. The mass func- 
tion drops quickly at very large masses, so that we can 
combine all masses above a certain limit to a single bin. 
We found it reasonable to define at any redshift 7 the cell 
boundaries by M n — 10™ h^ 1 M Q and the cell centers 
by M n = 10 n+1 / 2 /i _1 M 0) with M cut = M w . Because 
f(M > Mi 6 ,z) ~ to a very good approximation, it is 
sufficient to bin with an integer n between 10 and 15: 



M w 


< 


M 


< Mn 


"companion galaxies" 


M n 


< 


M 


< Mi 2 


"spiral galaxies" , 


M12 


< 


M 


< M13 


"elliptical galaxies" , 


Mis 


< 


M 


< M14 


"groups" , 


M x4 


< 


M 


< M X5 


"clusters" , 


AfiB 


< 


M 


< oo 


"superclusters" . 



The fraction of total mass in a bin is then 

rM n+ i 

A/„= / /dM, (72) 
Jm„ 

which correctly gives the total halo mass fraction (see 
Eq. (20)) when summed over all the halo bins: 

^A/ n =A/ H . (73) 

n 

Finally, the binned comoving halo densities are 

An n = ^Af n , (74) 

which also allows to relate the mass M n to the interhalo 
scale X n = Ann 1 ^ 3 - 



7 A more efficient binning may be obtained by adopting redshift- 
dependent mass bins M n (z). One could, for example, use bins 
such that AfLu has always the present-day value. See Section 
II B for a discussion about A/^jy. 



12 




FIG. 3: This illustration shows how we project the cylindrical 
filaments on the plane perpendicular to the geodesic. The 
triangular section on the right is moved to the left end of the 
filament so that the projected profile is independent of le- 




FIG. 4: Sketch of a cylindrical filament projected on a plane 
perpendicular to the geodesic. The effective length of the 
cylinder is Le of Eq. (77) as shown in Fig. 3. 



The halo impact parameter bh in the lens plane is re- 
stricted to the range bh € [0, R], where R is the comoving 
halo radius. It should be noted that through Eq. (33) the 
radius depends both on the halo mass and the redshift. 
See Fig. 2 for an illustrative sketch. It is efficient to dis- 
cretize b into bins of constant integrated surface density 
(constant equal mass): 

/ db2nbZ h (b,M n ,t t )= (75) 

Jb* nm N B 

where Ng is the number of 6-bins used and the surface 
density is given by (see Fig. 2): 

r R 2xdx 

X h {b h ,M,t) = M -=== <p h {x,t). (76) 
Jb h \A - K 

After the redshift and the halo mass is specified, 
the bin boundaries b\ nm and the weighted centers of 
gravity b^ nm can be computed using Eq. (75). Af- 
ter these are defined, the corresponding area functions 
AA in,n = T&inm+i - w &i n m, needed for the Poisson pa- 
rameters AN^ nm and the binned surface densities T% nm = 
£ h {b% nm ,M n ,ti) = M n /(N%AA!> nm ), can be computed. 

Overall, our halo model is described by the redshift in- 
dex and two internal indices for the discretized mass and 
impact parameter. In terms of the generic parameter u 
we can formally express this as: {h, inm} G {iu}, where 
the label h refers to the "halo" family of the inhomo- 
geneities. 



our filaments objects are confined to within a given Ar^ 
slice, they generically have only five degrees of freedom. 
Special symmetries of a sphere reduced the number of 
relevant parameters for halos to a sigle one, the impact 
parameter bh . For filaments the situation is slightly more 
complicated because a cylinder is invariant only for rota- 
tions around one symmetry axis. This leaves us with four 
degrees of freedom. The relevant parameters, however, 
are reduced by the symmetries to only two: the angle 9 of 
the filament main axes with respect to the geodesic and 
the impact parameter bf. For illustration see Figs. 3-4. 
This can be understood by the fact that all the possible 
filament configurations with the same 9 and bf have the 
same surface density and so can be resummed using 
the theorem given above. In particular, the projected 
surface density of our filaments does not depend on the 
coordinate lg along the projected main axes 8 allowing to 
rcsum the Zg-bins into a single effective bin of the total 
^-dependent length of the projection: 

L e = L sin 9. (77) 

The angle 9 is physically important because objects seen 
with a smaller 9 have a smaller cross section and a higher 
surface density. Accounting for this geometrical effect 
leads to a more strongly skewed lensing PDF. 



B. Discretizing the Filament Parameters 

The position and the orientation of a generic object 
in the three dimensional space is described by three po- 
sition and three orientation degrees of freedom. Since 



In order to have the projected profile independent of lg we 
"move" the triangular section as shown in Fig. 3. Strictly speak- 
ing our filaments do not have a cylindrical symmetry, but a 6- 
dependent "tilted" rotation symmetry, such that the projected 
surface densities follow equation (78) for all le- This approxima- 
tion should be good for long filaments with L>ii. 
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The expression for the surface density is given by density objects as: 9 

^ „ „ x M f f R 2xdx . . 



Ar / = — L_ y f 



(80) 



For the uniform filament profile this gives simply Y,f = 
2pfg s jR 2 — bj, where p/g = M f /(irR 2 Lg). For the gaus- 
sian profile (44) one instead finds 



/ i7T 

£, = 0.997-^ iJi/— exp[-^/2(i2/3) 2 ]. (79) 



The angle can be restricted by symmetry to the inter- 
val 9 £ [0, 7r/2], which we divide into At uniform length 
bins around centers 9 t = ir(t — \)/(2Nt)- The impact 
parameter is restricted to 6/ € [0,R], and it will be dis- 
cretized into Ng bins using the equal mass criteria, as in 
the halo case defined in Eq. (75). That is, our filaments 
will be represented by a family of "bars" of different sur- 
face densities T,{ tm and length L t = Lsin# t . The surface 
area of these bars to be used in the Poisson parameters 

Overall, our filament model is described by the redshift 
index, and two internal indices for the discretized angle 
and impact parameter. In terms of the generic parameter 
u we can again formally express this as: {/, itm} € {iu}, 
where the family label / refers to the "filament" family 
of the inhomogeneity. Of course adding, e.g., a mass 
distribution for filaments and/or an Z-depcndent density 
profile would increase the necessary number of indices in 
the discrete model. 



VI. CONFINING HALOS TO FILAMENTS 

Until now we have implicitly assumed that the halo 
and the filament families are independent. We will now 
show how a more realistic model with halos confined to 
the filaments can be set up in the sGL approach. The 
central observation is that the exact positions of halos 
and filaments within the co-moving distance bins are ir- 
relevant in Eq. (55) for the convergence. It then does not 
matter if we really confine the halos to filaments; for the 
same effect it is sufficient to merely confine them into the 
equivalent volume occupied by the filaments in a given 
bin. 

We can actually do even better than imposing a simple 
volume confinement. Indeed, it is natural to assume that 
the halo distribution follows the smooth density profile 
that defines the filament, which can be accounted for by 
weighting the volume elements by their reduced density 
profiles. Because only the projected matter density is 
relevant for lensing, this weighting is accounted for by 
the effective surface densities of the objects. We then 
introduce the effective co-moving thickness for the low 



where Pfi = Mf/Vfi is the average mass density of the 
filament with volume Vfi. The total effective length of a 
configuration of filaments {k( v } along the line of sight is 
then: 

Ar/({fc/J) = k l^{ v = An q fl . (81) 



It is easy to show that the statistical average of Ar( 
over the configuration space is just the expected distance 
covered by filaments: 



(82) 



where we defined the comoving average filament volume 
fraction: 



Qfi = Vfi An fi ■ 



(83) 



The confining simulation can now be performed in two 
steps: in the first step a random configuration of fila- 
ments is generated and Eq. (81) is used to compute the 
effective lengths Ar( at different co-moving slices. These 
lengths define then the Poisson parameters for the halos 
through Eq. (52): 



AA^({^})=A<Ar/A< 



(84) 



Note that the halo densities AnQ are corrected to ac- 
count for the reduced volume available due to confine- 
ment: 



An{ 



An h 



(85) 



where An^ u are the usual halo densities of Eq. (74) when 
the halos arc uniformly distributed through all space. In- 
serting Eq. (85) back to (84) and using Eqs. (81-82) we 
can rewrite Eq. (84) simply as 



ANf* 



Qfi 



AN h = X fi AN 



h 

iu 7 



(86) 



where AA^ arc the usual configuration independent 
Poisson parameters computed without confinement. 

Note that the central quantities needed in the evalu- 
ation of Eq. (86) are the filamentary weights AN- V and 



9 Note that this discussion is again general, and applies to all types 
of low density objects, not necessarily with cylindrical geometry. 
To distinguish filaments from halos we use v to refer to generic 
internal filamentary degrees of freedom. 
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surface densities E£ . With these given, all the infor- 
mation about the specific low contrast objects used is 
condensed into the single stochastic variable: 



Z-^V IV IV 

T AN f Y, f 

Z— 4V IV IV 



(87) 



which has an expectation value of unity and a mode 
smaller than unity. The lensing effects due to large-scale 
structures are thus tied to the probability distribution of 
Xf \ its added skewness is the effect coming from confin- 
ing halos within filaments. This opens the possibility to 
investigate and compare different filamentary geometries 
by means of the X/-PDF. We will develop these thoughts 
further in a forthcoming paper [32]. 

We can further generalize the picture given in this Sec- 
tion by considering different levels of confinement for ha- 
los of different masses. Indeed, because halos of different 
mass are treated independently in Eq. (86), one could 
have small halos populate also the voids and massive ha- 
los only the filamentary structure. 

Let us finally point out that in practical computations 
the convergence PDF is obtained by creating a large num- 
ber of independent halo configurations for each "master" 
filament configuration. Typically we use simulations with 
few hundred master configurations with a few hundred 
halo configurations each. 



addition to the lensing PDF. The power spectrum can 
indeed place useful constraints on the filament parame- 
ters which, as remarked in Section II B, are not tightly 
constrained by observations. 

While the power spectrum is relevant for understand- 
ing the correlations at the largest observable scales, the 
lensing PDF depends mainly on the much stronger in- 
homogeneities at smaller scales. This can be understood 
from Eq. (15) which shows the direct dependence of the 
lensing convergence on the density contrast. The small 
lensing variance from the large-scale correlations induced 
by the linear power spectrum was numerically computed, 
for example, in [41]. Moreover, weak lensing and the 
power spectrum probe somewhat different aspects of the 
inhomogeneities. The web-like structures of filaments 
and voids that affect weak lensing are mainly described 
by higher order correlation terms beyond the power spec- 
trum, and so special care has to be put in designing the 
filamentary structure. This is indeed the goal of the sGL 
method: to give an accurate and flexible modeling of the 
universe as far as its lensing properties are concerned. 



VII. RESULTS 



A. Power Spectrum 

In this Section we will discuss the power spectrum of 
our model universe and the importance of large-scale cor- 
relations for the lensing PDF. We start by making a con- 
nection between the sGL method and the so-called "halo 
model" (see, for example, [33-39] and [40] for a review), 
where inhomogeneities are approximated by a collection 
of different kinds of halos whose spatial distributions sat- 
isfies the linear power spectrum. The idea behind the 
halo model is that on small scales (large wavenumbers 
k) the statistics of matter correlations are dominated by 
the internal halo density profiles, while on large scales the 
halos are assumed to cluster according to linear theory. 
The two components are then combined together. 

The sGL modeling of the inhomogeneities can be 
thought as a two-step halo model where we first create 
the random filamentary structures and then place the 
halos randomly within these structures. Similarly to the 
halo model, one can then combine the linearly-evolved 
power spectrum with the nonlinear one coming from the 
filaments and the halos they contain. In this sense our 
filamentary structures extend the halo model by intro- 
ducing correlations in the intermediate scales between 
the halo substructures and the cosmological scale con- 
trolled by the linear power spectrum. In the halo model 
the power spectrum can be computed analytically, but 
here the calculation has to be done numerically. We 
are currently extending the sGL method such that the 
simulation will produce also the power spectrum [32] in 



Our main goal in this Section is to compare the sGL 
approach with the convergence PDF computed from large 
scale simulations. The idea is that by achieving a good 
agreement with numerical simulations we are proving 
that the sGL method does provide a good and accurate 
description of the weak lensing phenomena. As we shall 
see, this is indeed the case, as we can naturally reproduce 
the lensing PDF of the Millenium Simulation. 10 

Let us stress that while a comparison to simulations is 
a good benchmark test for the sGL approach, the simu- 
lations and underlying ACDM model do not necessarily 
describe Nature. Indeed observations do not yet provide 
strong constraints to the lensing PDF, which leaves room 
for very different types of large scale structures, examples 
of which have been studied for example in Rcf. [1, 10, 42]. 

With the accuracy of the method tested, one can re- 
liably compute the effects of selection biases using sGL. 
While we lack the necessary expertise to quantitatively 
estimate the possible bias parameters, we will study a 
simple toy model for the survival probability function 
P?" r to show qualitatively how selection effets might bias 
the observable PDF. Finally, we will also show how a 
JDEM-like survey could constrain the lensing PDF rel- 
ative to a given cosmological model. We shall begin, 
however, by the comparison with the simulations. 



We plan to compare the sGL model with other simulations in 
order to check its accuracy for different cosmological parameters. 
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Lensing PDF 

TABLE I: Filament parameterizations. 10 



Parameters I II 

A/o 4.5 9 

R p0 (h^Mpc) 2 4 

L p0 (ft _1 Mpc) 20 25 

profile uniform gaussian 

Mf (/i- 1 fi M oM ) 9.3 x 10 14 9.3 x 10 1E 



A. Comparison with the Millennium Simulation 

We shall now confront the sGL method with the 
cosmology described by the Millennium Simulation 
(MS) [14]. Accordingly, we will fix the cosmological pa- 
rameters to h = 0.73, Qjwo = 0.25, Hqq = 0.75, w = —1, 
Qbo = 0.045 erg = 0.9 and n s = 1. Moreover, the mass 
function of Eq. (31) agrees with the MS results. These 
parameters completely fix the background cosmology and 
the halos. 

We are then left with the filament parameters to spec- 
ify. First, we fix j3f = 0.5, which means that half of 
the unvirialized mass forms filaments, while the other 
half is uniformly distributed. This value determines the 
maximum possible demagnification Amu e in the lensing 
PDF, corresponding to light that misses all the inhomo- 
geneities. For the background model described above this 
implies Attiue — 0.17 mag. 

The next parameter to set is the present-day filament 
mass Mf - Since j3f was taken a constant, a lower Af/rj 
implies a higher comoving filament density (see Eq. (41)), 
a more homogeneous universe and hence less pronounced 
lensing effects. In order to connect with the literature 
we actually have to use the dressed filament mass of 
Eq. (45) which, because of the chosen /3f, is roughly 
Mj? w 3-M/o- To choose the parameters defining Mfo, 
we follow Ref. [43], according to which it is reasonable 
to fix the average filament overdensity to Aw = 4.5 
and the filament radius to R p o = 2h~ 1 Mpc. Note 
that the overdensity of the dressed filament is roughly 
3A/o and that our filaments include also the large ha- 
los at the filament ends. Moreover the density profile 
seems rather uniform within R p o according to [43] and 
so we use the simple uniform density profile. The dis- 
tribution of filament lengths is quite wide (see Fig. 7 of 
Ref. [43]) and we choose to use a representative value of 
L p o = 20/i~ x Mpc. With these choices the bare filament 
mass becomes Mfo = 3.1 x 10 14 ]i~ 1 CImoMq. We have 
summarized in Table I this first filament parameteriza- 
tion (I) . We can now compare the properties of the lens- 
ing PDF from the sGL method to the one from the Mil- 
lenium Simulation, which is a courtesy of Refs. [11, 12]. 

The four histograms in Fig. 5 show the four steps in 
building the model universe in the sGL method. The 
first (orange) histogram from the top corresponds to a 
universe endowed only with the halos specified by the 
mass function. The second (green) histogram from the 
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FIG. 5: Shown are the lensing PDFs for a source at z — 1.5 
for the ACDM model of the Millennium Simulation [14]. The 
dashed line is the lensing PDF generated by shooting rays 
through the MS [11, 12], while the histograms are obtained 
with the sGL method. The first (orange) histogram from 
the top corresponds to a model universe with halos but no 
filaments, the second (green) is relative to a universe with ha- 
los and cylindrical filaments randomly placed and the third 
(blue) to the case in which the halos are confined within mass- 
less filaments. The fourth (red) histogram is the most realis- 
tic case with halos confined within randomly placed massive 
filaments. The filaments are modelled according to Table I 
(parameterization I). See Section VII A for details. 



top-left represents a model in which half of the uniform 
matter density (/?/ = 0.5) is condensed into filaments, 
but halos are not confined to them. The third (blue) 
histogram from the top-left represents a model in which 
the halos are confined within the filaments but these are 
massless. Finally, the lowest (red) histogram represents 
the case where the halos are confined within massive fil- 
aments. 

For these parameters the lensing PDF seems domi- 
nated by the halo contribution due to the NFW density 
profiles and the effect of the large scale clustering seems 
to increase the dispersion. In particular, the clustering of 
halos and of smooth mass within filaments seems to give a 
similar contribution as the two corresponding histograms 
are close to each other, with a slightly more skewed PDF 
for the former as compared to the latter. In other words, 
the nonlinear clustering of halos at the filament scales 
appears to increase the dispersion similarly to the distri- 
bution of smooth matter in low density filamentary struc- 
tures. The results relative to the Millennium Simulation 
from Ref. [11, 12] are shown with a dashed line. The 
agreement becomes very good when the large scale struc- 
tures are fully included, showing that the lensing PDF is 
a sensitive probe of the large scale structures. We can 
conclude that the halos have to be confined within mas- 
sive filaments if we are to reproduce the MS results, at 
least for the parameters chosen. Moreover, by fixing the 
other parameters one could place an upper bound on fila- 
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case. The poor agreement, however, does not mean that 
the universe cannot contain such large structures, but 
merely that the MS does not have them. We will further 
discuss this example in the next Section. 

We have also compared the redshift dependence of the 
mode and the dispersion of the lensing PDF with the 
MS results [11, 12] using the filament parameterization 
I. We choose the full width at half maximum (FWHM) 
as an indicator for the dispersion instead of the standard 
deviation a and, to be precise, we will compute the quan- 
tity v = FWHM/2.35. The motivation for this is that v 
is insensitive to the long high-magnification tail and so 
represents a more robust estimate of the dispersion in 
the region where the lensing PDF is effectively nonzero 
(i.e., in the very weak lensing regime). Moreover, for a 
gaussian PDF one finds v = a. The results are shown in 
Fig. 6, and the agreement with the results of Ref. [11, 12] 
is again very good. These results are very encouraging 
because we next want to study quantitatively the selec- 
tion effects using the sGL method. 



FIG. 6: Shown is the redshift dependence of mode and dis- 
persion v = FWHM/2.35 of the lensing PDF for the ACDM 
model of the Millennium Simulation [14]. The datapoints 
were obtained with the sGL method (parameterization I of 
Table I), while the dashed line refers to lensing PDFs gen- 
erated by shooting rays through the Millennium Simulation 
(MS) [11, 12]. 



ment masses. Indeed, by considering the same modeling 
of parameterization I but with j3f = 1 and AfQ = 9, we 
have a universe model with the same geometrical prop- 
erties 11 but with twice as massive filaments. The result- 
ing lensing PDF has now a longer low-magnification tail 
(Am > 0) and is not in good agreement with the MS 
result. Similarly, we could look for degeneracies between 
background parameters and filament parameters. For ex- 
ample, still considering the parameterization I, we have 
found that a lower value of erg = 0.8 is compensated to a 
good degree by using a gaussian filament profile instead 
of a uniform one. 

The sGL method is an ideal tool to explore the changes 
induced in the lensing PDF for different choices for the 
filament parameters and in Fig. 7 we show the lensing 
PDF for a universe in which the filaments are ten times 
more massive. The same parameters were used as before 
but now with A/ = 9, R p q = 4/i _1 Mpc and L p0 = 
25fo~ 1 Mpc. Moreover the filaments have the gaussian 
profile of Eq. (44) and (79). These parameters are again 
summarized in Table I (parameterization II). It is clear 
that the agreement with the MS results is lost in this 



11 The average filament volume fraction is unaltered if the ratio 
(3f/Afo is kept constant. 



B. Selection Bias Effects 

We now introduce a toy model for the selection effects 
to evaluate qualitatively their impact on the lensing PDF. 
We model the survival probability simply as a step func- 
tion in the impact parameter: 

P sur (6) = {? fe /f <Scut , (88) 
I 1 otherwise 

where s cu t gives the opaque fraction of the halo radius. 
In Fig. 7 we show the lensing PDF for the filament pa- 
rameterization II in the case of s cut =0.1 (dotted PDF). 
The lensing PDF with unitary survival probability is also 
shown for comparison (solid PDF). As expected the se- 
lection effects reduce the high magnification tail of the 
PDF. The effect is more clearly seen by plotting the dis- 
persion as function of redshift for the two cases. The 
top panel of Fig. 8 shows indeed that the standard devi- 
ation cr is appreciably smaller when the selection effects 
are turned on. The dispersion v (middle panel) and the 
mode (bottom panel) are, instead, much less affected. 

The conclusions are twofold. First, this shows that 
v is indeed more robust than a with respect to the high 
magnification tail, which is effectively cut by the selection 
effects. Moreover, the value of a seems to move to the 
value of v when the survival probability is nontrivial. 

Second, while selection effects may reduce the variance, 
larger scale inhomogeneities have the opposite tendency 
and the two can compensate each other. We illustrate 
this issue by also plotting in Fig. 8 the results relative to 
parameterization I without selection effects. Clearly the 
more massive filaments of the parameterization II with 
selection effects give an effective a very close to the one 
of the lighter filaments of the parameterization I without 
selection effects. In other words, for a given observational 
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FIG. 7: Lensing PDF for a source at z = 1.5 for the filament 
parameterization II of Table I. The solid (green) histogram 
is without selection effects, while the dotted (blue) histogram 
is with the selection effects toy-modelled by Eq. (88). For 
comparison with Fig. 5, the lensing PDF relative to the Mil- 
lennium Simulation is plotted again as a dashed line. See 
Section VII for details. 



bound on the lensing variance, there may be some degen- 
eracy between observational biases and weak lensing by 
large scale structures. This clearly shows the importance 
of correctly modeling selection effects. The above degen- 
eracy, however, is broken by the redshift dependence of 
v and mode, showing that a precise measurement of the 
lensing PDF (beyond the variance) can yield important 
cosmological information. We will discuss in the next 
Section the possibility of observing the lensing PDF. 
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FIG. 8: Shown is the redshift dependence of a (top panel), 
v = FWHM/2.35 (middle panel) and mode (bottom panel) for 
the filament parameterization II of Table I with (dotted line) 
and without (dashed line) the selection effects of Eq. (88). 
The solid line is relative to the filament parameterization I 
without selection effects. 



C. Measuring the Lensing PDF 

We will now discuss the possibility of measuring the 
lensing PDF for the universe of the Millennium Simu- 
lation. We will perform a simplified analysis in which 
we assume that we are observing perfect standard can- 
dles and we are thus neglecting the intrinsic dispersion 
in the SNe absolute luminosity. Moreover, we will also 
neglect any observational uncertainties. We refer to sub- 
sequent work for a more comprehensive analysis, while 
here we give an illustrative picture for the prospects of 
observationally constraining the lensing PDF. We will in 
particular consider a dataset of a JDEM-like survey with 
one thousand high redshift SNe. With this in mind, we 
show in Figs. 9-10 the lensing PDFs for sets of 200 and 
1000 perfect standard candles at a given target redshift 
of z = 1.5. Also plotted are the 2a errors, where the stan- 
dard deviation was computed with turboGL in a Monte 
Carlo fashion by generating many of such histograms and 
calculating the a for each bin height. 

The two figures represent two possible comparisons of 
the data with a given cosmological model, the difference 



being the amount of prior cosmological information as- 
sumed. Indeed, since the observational data is heteroge- 
nous it must be binned to some redshift intervals Az 
and the SNe luminosities have to be normalized to the 
same bin redshift. If Az is small the implicit dependence 
of the above correction on the background cosmology is 
small and one finds a direct constraint on the cosmo- 
logical model. This situation is displayed in Fig. 9 for 
a JDEM-type dataset with 200 SNe in the redshift bin 
Az = [1.45, 1.55]. It is clear that this single redshift bin 
cannot put very tight constraints on a generic lensing 
PDF and one has to use all the high-redshift dataset. 

One way to proceed is to compare per bin the predicted 
lensing PDF with the observed one, that is, by treating 
different bins independently. Another way could be to 
combine all the 1000 SNe at z > 1 to one bin. Fig. 10 dis- 
plays the relative histogram, which now has much smaller 
error bars as compared to Fig. 9. The normalization of 
the dataset to the same target redshift, however, used 
heavily the prior assumptions about the cosmological 
model. That is, the comparison shown in Fig. 10 pro- 
vides only a consistency check for the assumed cosmolog- 
ical model, and it should be noted that the reconstructed 
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FIG. 9: Lensing PDF for a source at z — 1.5 for the filament 
parameterization I of Table I (without selection effects). The 
histogram simulates a possible PDF obtained with a dataset 
of 200 perfect standard candles. The 2a errors are represented 
as rectangles. The PDF relative to the Millennium Simulation 
is plotted as a dashed line. See Section VII C for details. 



Lensing PDF 




FIG. 10: Same as in Fig. 9 but for a dataset of 1000 perfect 
standard candles. 



lensing PDF might hide a high degree of degeneracy on 
the cosmological parameters. 



VIII. CONCLUSIONS 

In this paper we have revised and extended the stochas- 
tic approach to cumulative weak lensing first introduced 
in Ref. [1]. The (sGL) method is based on the weak 
lensing approximation and on generating stochastic con- 
figurations of halos and large-scale structures along the 
line of sight to cosmological objects, such as supernovae. 
The halos are intended to model the virialized structures, 
such as galaxies and clusters of galaxies. We model these 



structures accurately by means of realistic mass distri- 
bution functions and by NFW density profiles (with ap- 
propriate z and M dependent concentration parameter). 
The filamentary structures were modelled with (possibly) 
non-uniform cylindrical objects. 

This modeling incorporates the essential fact that 
most of the spatial volume is occupied by underdensitics 
(voids), while most of the mass lies in overdense regions of 
much smaller total volume, in the form of virialized clus- 
ters and large filamentary structures. As a consequence 
the column mass density along a random single geodesic 
is likely to be lower than the average, which gives rise to 
the skewness of the lensing PDF. A skewed PDF is of po- 
tential importance in particular for the interpretation of 
the supernova observations, which are still probing much 
smaller angular scales than the scale at which the homo- 
geneity is recovered. The sGL method with its excellent 
performance provides a useful tool to investigate the im- 
pact of lensing, especially on analyses based on likelihood 
approaches as done in Ref. [10]. 

We confronted our sGL modeling with the lensing re- 
sults [11, 12] from the Millennium Simulation [14]. Even 
though the MS does not necessarily represent the actual 
universe, this comparison provides a very useful bench- 
mark test, because here the results are exactly known, as 
well as the main characteristics of the inhomogeneities 
involved. The sGL-generated PDF agreed with the MS 
results very well and naturally without any fine tuning, 
when all model parameters were chosen to match the 
typical structures seen in numerical simulations. In the 
MS cosmology most of the lensing was found to be due 
to the halos, with a smaller, but clearly distinguishable 
contribution coming from filaments. We conclude that it 
is necessary to embed the halos along massive filaments 
to reproduce the MS results. 

Moreover, the sGL method can easily incorporate ob- 
servational biases such as selection effects or any other 
mechanisms that can obscure lines of sight. We believe 
this to be an appealing feature of the sGL method, es- 
pecially with the growing size of SNe datasets: a large 
number of observations indeed reduces the statistical bias 
due to the skewed lensing PDF, but does not eliminate 
the systematic ones possibly caused by selection effects. 
One can alleviate such problems with a good understand- 
ing and modeling of the biases, which is also needed if 
one wants to observationally constrain the lensing PDF. 
Indeed, we have shown by a simple toy-model example 
that selection effects can potentially alter the observable 
PDF in such a way as to cancel the opposing modification 
coming from very large scale structures. Assuming no se- 
lection effects might then lead to a premature exclusion 
of such structures. 

Finally, we have also shown how a JDEM-like survey 
could constrain the lensing PDF relative to the cosmol- 
ogy of the Millennium Simulation. Our results suggest 
that such a survey might be able to give important in- 
formation about the PDF and thus about the large-scale 
structure of the universe. 
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Along with this paper, we release an updated version 
of the turboGL package, a simple and very fast Mathc- 
matica implementation of the sGL method, which can be 
found at turboGL.org. 
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Appendix A: Performance of the sGL Method 

In order to estimate the performance of the sGL 
method we have to calculate the total number Nb ins 
of {m}-bins. It is easy to find that Nuns/NsNs — 
Nh + NpNr, where Nh is the number of halo-mass 
bins, N F is the number of filaments, N T is the num- 
ber of 0-bins and we have used for simplicity the same 
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number of impact-parameter bins Nb for both halos and 
filaments. We remind that Nt is the effective number 
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